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Triply deuterated ammonia in NGC 1333 
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Abstract. The Caltech Submillimeter Observatory has detected triply deuterated ammonia, ND3, through its 
Jk = 1§ — *• 0o transition near 310 GHz. Emission is found in the NGC 1333 region, both towards IRAS 4A and a 
position to the South-East where DCO + peaks. In both cases, the hyperfine ratio indicates that the emission is 
optically thin. Column densities of ND3 are 3 — 6 x 10 11 cm -2 for T cx =10 K and twice as high for T ex =5 K. Using 
a Monte Carlo radiative transfer code and a model of the structure of the IRAS source with temperature and 
density gradients, the estimated ND3 abundance is 3.2 x 10~ 12 if ND3/H2 is constant throughout the envelope. In 
the more likely case that ND3/H2D + is constant, ND3/H2 peaks in the cold outer parts of the source at a value 
of 1.0 x 10 -11 . To reproduce the observed NH3/ND3 abundance ratio of ~1000, grain surface chemistry requires 
an atomic D/H ratio of «0.15 in the gas phase, >10 times higher than in recent chemical models. More likely, 
the deuteration of NH3 occurs by ion-molecule reactions in the gas phase, in which case the data indicate that 
deuteron transfer reactions are much faster than proton transfers. 
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i-h ! 1. Introduction 



Deuterium-bearing molecules have attracted attention in 
recent years. Physically, they appear to be good probes of 



H and CO onto grains, followed by reaction, is thought 
to produce solid H2CO and CH3OH. This mechanism 
favours deuteration because the atomic D/H ratio in the 
gas phase is much greater than the elemental ratio. The 
observed abundances of HDCO, CH3OD and D2CO in 
Orion and IRAS 16293 indicate their synthesis on dust 



the ven y cold phases of molecular clouds prior to star for- 
mation . Chemically, the isotopic composition of molecules 
is an important clue to their formation mechanism. There 
are tw p ways to make deuterated molecules. First, at 
temper atures ^70 K, the gas-phase reaction equilibrium 



grain s (Turner 199C; Charnley et al. 1997; Ceccarelli et al 



1998 ). The same mechanism may work for NH3, provided 



most nitrogen is in atomic form, and support for this route 
comes from detections of solid NH3 ( |Lacy et al. 1998 ; |Gibb 
et al. 2000). However, these detections remain tentative 



HD H2D" 1 " + H2 is shifted in the forward direc- 
tion. Subsequent deuteron transfer from H2D + to, e.g., 
CO and N2, leads to the large observed abundance ratios 
of DC O+ /HCO+ and N9D+/N9H+ of ~0.1. four orders 



( Dartois fc d'Hcndccourt 2001 ) , and observations of N2H" 1 " 
suggest that in dense clouds, most nitrogen is in molecular 
form ( |Womack et al. 1992j ). 



While NH2D has been observed in many sources (Saitc 



of magnitude higher than the Galactic D /H ratio (jTurncr| ct al - 2 000|; |Shah fc Wootten 200 1| ), ND 2 H has only been 



2001 ). The key species of this chemical scheme, H2D" 1 ", was 
r ecently detected i n the Class source NGC 1333 IRAS 4A 
( Stark et al. 1999 ). The high densities and low tempera- 
tures in this object promote the formation of H2D + out 
of H^ and HD, and also prevent its destruction because 
the major destroyer of H2D + , CO, is depleted by a factor 
of ^100 due to freeze-out onto dust grains. 

Alternatively, deuterium-bearing molecules can be 
formed on dust grains by surface chemistry. Accretion of 



detected so far towards the cold, starless cores L134N and 
L1689N QRoucff ct al. 2000|; |Loinard et al. 200 1|). The 



temperatures in these sources of «10 K are too low for 
significant evaporation of even the most volatile ices to 
occur. Gas-phase reactions can probably account for the 



observed abundance ratios (Rodgers & Charnley 2001) 
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but observations of ND3 would present a strong test, as 
the gas-phase route produces ~3 times more ND3 than 
the grain surface route. As a step towards measuring the 
relative importance of gas-phase and solid-state deutera- 
tion, we have observed the Jk — lo —* 0q transition of 
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Fig. 1. Top: Spectrum of NGC 1333 IRAS 4A near 309.9 
GHz, taken with the CSO, with the 3-component fit de- 
scribed in the text overplotted and the expected positions 
of the hyperfine components of the ND3 Jk — 1q — > line 
indicated. Bottom: Spectrum taken at the (+23, —06)" off- 
set position, where DCO + peaks. 

ND3 towards NGC 1333. To gether with sim ilar observa- 
tions towards Barnard 1 by Lis et al. (2002 ), this is the 
first detection of a triply deuterated molecule in interstel- 
lar space. 



2. Observations and Results 

The rotational energy levels of symmetric top molecules 
are labeled by the total angular momentum J and its pro- 



of 25", slightly larger than the diffraction limit. Initial 
observations of ND 3 J K = lg -> Of, toward NGC 1333 
IRAS 4A (a 200 o = 03 ,l 29 m 10 s .3, S 20 qo = +31°13'32") 
were carried out December 5-8, 2001. Additional data were 
taken January 24, 2002, with the local oscillator frequency 
shifted to verify the sideband origin of detected features. 
Weather conditions were average with zenith opacities of 
wO.l at 225 GHz. An IF frequency of 1594 MHz was used 
to put the Jk = 1q — * 0q transition in the image sideband. 
To subtract atmospheric and instrumental background, 
positions 240" offset were observed using the chopping sec- 
ondary mirror. The main beam efficiency at the time of ob- 
servations, measured through observations of Saturn, was 
64%, using a planetary brightness temperature of 135 K 
at 310 GHz. 

Figure |l| (top) shows the combined data from all 
nights, aligned in signal frequency. This spectrum repre- 
sents 151 minutes of on-source integration with a mean 
T sys of 521 K and a source elevation of 44°. A feature is 
detected at 309909.4 MHz, at which frequency neither the 
JPL (|Pickctt ct al 199Sj , spec.jpl.nasa.gov) nor the CDMS 
(Miillcr ct al. 2001, www.cdms.de) catalog lists any plau- 
sible molecular lines other than ND3. The 307 GHz com- 
ponent is not detected to T m b< 14 mK (la). 

Data were also obtained at a position 23" East and 
6" South of the IRAS 4A position, which is where the 
DCO + 3^2 emission peaks (Lis et al., in preparation). 
The bottom panel of Figure |l| shows that the ND3 line 
is also somewhat stronger than at the IRAS 4A position. 
The noise level for this spectrum is 23 mK (la). 

Lines of ND3 exhibit hyperfine structure due to cou- 
pling of the 14 N nuclear spin with the rotational angu- 
lar momentum. Splitting due to D, owing to its small 
quadrupole moment, is only «200 kHz and remains un- 
resolved in our data. Thus, we have fitted the observed 
spectrum with three Gaussian profiles with optical depth 
ratios of 5:3:1 and assuming that the components have 
equal widths and excitation temperatures. Such a fit has 
four parameters: optical depth r, central velocity Vb, line 
width AV and intensity T m b=r(T ex — Tb g ). We used the 
HFS method inside the CLASS package. Leaving all pa- 
rameters free gives the values reported in columns 2-5 of 
Table |l|. For the IRAS 4A position, fixing V and/or AV to 
values measured in other lines, 7.0 and 1.2 km s _1 (Blake 



jection on the molecular symmetry axis K. For NH3 and 
ND3, inversion motion splits each level further into states 
which are symmetric (s) and antisymmetric (a) upon re- 
flection in the plane of the H or D atoms. The measured 



et al. 1995) gives similar results. The fit results indicate 



I " 



On 



frequencies of the ND3 Jk 
309908.46 (F=l->1), 309909.69 (F=2 
(F=Q->1) MHz; those of the J K = lg 
306735.58 (F =! ->!), 306736.96 (F=2 
(F=Q-+1) MHz ( |Helminger ~ 
statistics, the Jk = 1q — > 0q transition is stronger than 
the Jk = lg — > 0q transition by a factor of 10. 

The Caltech Submillimctcr Observatory (CSO) is a 
10.4-m single-dish antenna located atop Mauna Kea, 
Hawaii. At 310 GHz, the CSO has an FWHM beam size 



transition are 
l") and 309911.53 
Og transition are 
1) and 306738.95 
Gordy 196E). Due to spin 



that the line shape is consistent with the optically thin 
hyperfine intensity ratio of 5:3:1, but that the signal to 
noise is not high enough to constrain the optical depth. 
Without taking hyperfine broadening into account, AV at 
the IRAS 4A position would be 2.6 km s _1 , much broader 
than other lines from the cold component of this source 
( Blake et al. 1995j ). The line shape thus confirms the as- 
signment of the line to ND3. 

From the observed line strengths integrated between 
3.5 and 10.0 km s" 1 , given in column 6 of Table [ij we esti- 
mate the ND3 column density assuming an excitation tem- 
perature T ex =10 K. In the optically thin limit, but outside 
the Rayleigh- Jeans regime, and using a background tern- 
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Table 1. Line parameters' 



Position 


T 


V 
km s _1 


AV 

km s _1 


mK 


J'T mb dV 
mK km s _1 


iV(T cx =10 K) 
10 11 cm" 2 


7V(T CX =5 K) 
10 11 cm" 2 


IRAS 4A 
DCO+ 


0.6(20) 
0.1(27) 


7.0(2) 
6.5(1) 


1.6(6) 
1.0(2) 


81(62) 
149(22) 


71(21) 
144(43) 


2.9(9) 
5.9(18) 


6.3(19) 
12.8(38) 



a Numbers in brackets denote the uncertainty in units of the last decimal 



perature of Tbg = 2.7 K, the velocity-integrated optical 
depth follows from 



T mh dV 
so that 

with 



hv ( 1 



1 



- 1 



dV 



N„ 



8irv 3 



AN u (e 



Jiv/kT h . 



hv/kT a , 



TdV 



1) 



9uN c - E ^/ kT ^ 



Q{T CX ) 



yields the column density estimates in column 7 
of Table pj Here, Q(T CX ) is the partition function 



-Ei/ 



equal to 38.3 for T cx =10 K. The 



Einstein A coefficient of 2.57 x 10 



follows from the 



dipole moment of 1.49 D (di Lonardo & Trombctti 1981). 

The assumed value of T ex represents a kinetic temper- 
ature at which chemical fractionation should be efficient. 
However, at densities well below the critical density of 
this line, 10 7 - 10 8 cm , T cx will drop below Tkin, which 
changes the column density estimate. As an example, the 
last column of Table [l] gives the values for T cx =5 K. 



3. Abundance of NDj 



To estimate the 
used the Monte 
gram by 



abundance of ND3 we have 
Carlo radiative transfer pro- 



Hogcrhcijde fc van der Tak (2000 



talisker.as.arizona.edu/^michiel/ratran.html). Lacking 
auxiliary data on the DCO + position, we concentrate on 
NGC 1333 IRAS 4A, for which we take the temperature 
and density structure from stark et al. (1999). Between 
the outer and inner radii of 3100 and 10 AU, temperatures 
increase from 13 to 320 K, and densities from 2 x 10 6 
to 4 x 10 11 cm" 3 ; AT(H 2 )=3.1 x 10 23 cm" 2 in a 13" 
beam, but strongly depends on beam size due to the 
R~ 2 density distribution. The radiative transfer model 
for ND3 includes the 30 terms up to 100 cm -1 above 
ground, including the inversion splitting but not the 
hyperfine structure. Rate coefficients for de-excitation 



of NH3 in collisions with H2 from Danby et al. (1988) 
arc used, scaled to the different reduced mass of the 
ND3-H2 system, and augmented with the terms that are 
Pauli-forbidden in NH3, and with transitions that would 
be ortho-para conversions in NH3. Initially, a constant 
abundance of ND3 (relative to H2) was assumed. The 



excitation of ND3 as a function of radius is calculated 
with the Monte Carlo program. The result is integrated 
over the line of sight and convolved with a 25" beam. The 
area under the synthetic line profile matches the observed 
value for ND 3 /H 2 = 3.2 x 10~ 12 . 

As an alternative model, the ND3 abundance was as- 
sumed to follow that of H2D" 1 ". The major chemical for- 
mation path to ND3 starts with the reaction of NH3 with 
H2D + and its derivatives DCO + and N 2 D + , and proceeds 
through NH 2 D and ND 2 H. As a simple way to model this 
behaviour, we have assumed a constant H2D + /ND3 ratio. 
However, this ratio would vary in the case of a varying 
NH3 abundance, and if the alternative route starting with 
the reaction of N + with HD competes, which is slightly 
endothermic. In the absence of sufficient constraints, we 
keep H 2 D + /ND3 constant. Our two assumed ND3 abun- 
dance profiles could be tested indirectly by observations 
of key deuterated molecules such as DCO + and N 2 D + . 

The H 2 D+ abundance profil e in NGC 1333 IRA S 4 A 
was calculated analytically by Stark et al. (1999 ), us- 
ing assumed values for the cosmic-ray ionization rate 
(C = 5xl0~ 17 s^ 1 ) and the abundances of HD (2.8 xl0~ 5 ) 
and D (2.8 x 10 -6 ), and using a CO abundance of 4 x 10~ 6 
estimated from C 17 data. Due to the small energy dif- 
ference between H3 and H 2 D + , the H 2 D + abundance is 
strongly peaked toward large radii where temperatures are 
low. We have re-calculated the H 2 D + abundance profile 
using £ = 2.6 x 10~ 17 s^ 1 , the mean of the values implied 
by observations of H^ and H 13 CO + towards seven mas- 



sive young stars (van der Tak & van Dishoeck 2000). This 



calculation also includes dust radiation which the one by 
|Stark et al. (1999| ) did not. The new calculations are still 
consistent with the measured H 2 D + line flux, and indicate 
an H 2 D + abundance increasing from 5 x 10~ 19 at a radius 
of 10 AU to 5 x 10~ 10 at R = 3100 AU. To model the ND 3 
data, models were run for several values of the H 2 D + :ND3 
ratio, and agreement between observed and calculated line 
flux was found for H 2 D+/ND 3 =46. The ND 3 abundance 
at large radii is then 1.0 x 10 -11 , a factor of 3 higher than 
that found assuming a constant ND3 abundance. 

4. Chemistry of ND 3 

Table |^ summarizes the measured column densities of NH3 
isotopomers toward NGC 1333 IRAS 4A. It is seen that 
Af(NH 3 )/Ar(NH 2 D) «10 and iV(NH 3 )/iV(ND3) «1000; no 
observations of ND2H exist yet. The available data sug- 
gest a trend where with each H— >D substitution, the col- 
umn density drops by an order of magnitude. Current 
models of gas-phase chemistry, on the other hand, predict 
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Table 2. Column densities of ammonia isotopomers to- 
wards NGC 1333 IRAS 4A 



Species 



N 



Beam Reference 



~NH^ 3.1(4) x 10 14 
NH 2 D 2.2(7) x 10 13 
ND 3 2.9(9) x 10 11 



74 
90 
25 



5hah fc Wootten (2001 ) 
5hah fc Woottcn (200 1| ) 
this work 



that N D 2 H/ND 3 > NH 2 D/ND 2 H > NH 3 /NH 2 D ([Rodgers 



& Cha rnley 2001 ), unless deuteron transfer reactions are 
much more rapid than proton transfers. The same trend is 
expected in the case of grain surface chemistry. However, 
the measured column densities are uncertain by «30% due 
to calibration, so their ratios could be off by a factor of 
two and cannot be used to rule out either mechanism. 

Since a straightforward comparison of column densities 
may be complicated by the differences in beam size of the 
data in Table ||, we have determined the NH3 abundance 
toward NGC 1333 IRAS 4A using the approach of § |. 
The same temperature and density structure as in § H arc 



used, and the original collisional rate coefficients of Danby 
et al. [198S|). Based on the UMIST database flLc Tcufi 



et al. 2]00C| , www.rate99.co.uk), the rates of the major de- 
struction reactions of NH 3 do not depend on temperature. 
Formation of NH 3 is mainly by dissociative recombination 
of NH4 , the rate of which has a T~ - 5 dependence, which 
in the model for NGC 1333 IRAS 4A corresponds to a fac- 
tor of 5. This factor is dwarfed by the exponential increase 
in H 2 D + so only constant-abundance models have been 
considered for NH 3 . The observations of [3hah fe Wootter] 
(200lfi]are reproduced for NH 3 /H 2 = 1 x 10~ 8 . 

Rodgers & Charnley (2001) present a chemical scheme 
to form deuterated ammonia in the gas phase, which 
assumes that the branching ratios of dissociative re- 
combination are statistical and that the relevant reac- 
tion rates are isotope-independent. Using this scheme, 
for NH 3 /NH 2 D=10 (Table §), an abundance ratio of 
NH 3 /ND 3 of wlOOOO is expected. The observed value of 
10 _8 /10 -11 = 1000 is inconsistent with this prediction. 
This disagreement may indicate that a detailed chemical 
network is needed instead of a statistical treatment. In ad- 
dition, NH 3 /NH 2 D was measured in an arcmin-sized re- 
gion and may be < 10 within the 25" CSO beam. Based on 
the observed NH 3 /ND 3 ratio, NH 3 /NH 2 D could approach 
unity on small scales. 

In the case of surface chemistry, the observed 
NH 3 /ND 3 ratio implies an atomic D/H ratio of «0.15 
in the gas phase (Elodgers & Charnley 2001). This is 
significantly higher than the values of 10~ 2 — 10~ 3 in 
the chemical models of Roberts fe Millar (2000| ). Based 
on this discrepancy and on the closer agreement of the 
observed NH 3 /ND 3 ratio with the gas-phase prediction 
for NH 3 /NH 2 D=10, we tentatively conclude that ion- 
molecule reactions are presently the preferred formation 
mechanism of NH 3 in NGC 1333. This mechanism can re- 
produce the observed deuteration levels if deuteron trans- 



fer reactions are much faster than proton transfers. In the 
future, this conclusion should be tested through measure- 
ments of the NH 3 , NH 2 D, ND 2 H and ND 3 abundances in 
a larger source sample. Any conclusion drawn from such 
data will be much stronger if the lines are measured with 
similar beam sizes. More detailed chemical networks are 
also needed. Such a project could constrain the relative 
importance of gas-phase and grain-surface deuteration as 
a function of environment. 
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